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Abstract. Properties of transactinide nuclei needed for the interpretation of spectroscopic 
data and for the extrapolation to the expected "island of stability" of superheavy nuclei are 
reviewed from a theorist's point of view in the framework of self-consistent mean-field models. 



1. Introduction 

Superheavy elements (SHE) are usually defined as the chemical elements whose atomic nuclei 
have a vanishing fission barrier in the liquid drop model. In this schematic picture, for nuclei with 
Z > 106 the strong Coulomb repulsion pushes the nucleus to large deformation and overcomes 
the surface tension that tries to keep the nucleus at small deformation. That such systems can 
exist and have been indeed observed is due to quantal shell effects: according to Strutinski's 
theorem, a less than average density of single-particle levels around the Fermi energy gives 
additional binding, whereas a larger than average value reduces the binding. As a consequence, 
the variation of the single-particle spectrum with deformation creates local pockets and barriers 
in the deformation energy surface. Being on a scale of about 10 MeV, the amplitude of these 
variations is small compared to the total binding energy of almost 2 GeV, but it is large enough 
to stabilise many heavy nuclei against spontaneous fission, a-decay, /3-decay or electron capture 
often become their dominating decay modes instead. SHE are also often called tmnsactinides. 

Over the last four decades, the limits of the chart of nuclei have been pushed far into the 
realm of superheavy elements up to Z = 118 HI El |3l IH [5l El [7]. One motivation for the 
experimental efforts in this field of nuclear research is the long-standing prediction of an island 
of very long-living spherical SHE, which originally was expected to be centred at .Z^ ~ 114 and 

~ 184 [HI [9l [TOl [TT| I12j . Even when the expectations about their lifetimes had to be cut 
down from years or even longer to minutes, perhaps just seconds, finding such nuclei would still 
provide a unique laboratory to study the nuclear many-body system under extreme conditions. 

Because of the quickly dropping fusion cross sections [H EJ [3l HI [5l [13] , the number of events 
for the production of elements with Z > 104 is very low, ranging from just a few to a few dozen. 
Consequently, the information available for these nuclei is limited to their existence (where even 
the attribution of a neutron number is sometimes a matter of debate), their dominating decay 
channel (usually fission or a decay) and corresponding lifetime, and in case of a-decay also the 



Qa value. For some odd-74 nuclei the branching of a-decay into different levels in the daughter 
nucleus is observed already at this low level of statistics. Not only the properties of the atomic 
nuclei are studied for these heavy elements, but also some of their basic chemical properties |14j . 

Detailed spectroscopic information about their nuclear structure is available only for nuclei up 
to the lower border of the region of transactinides, for Z values up to about 104 [TS t fTG t fTT l fTS t fTQ]. 
These can be produced in much larger quantities and have longer lifetimes than most of the 
(known) heavier nuclides, such that they can be studied with a wide range of experimental setups. 
Electromagnetic transitions within and between rotational bands up to quite high spin are 
obtained from in-beam spectroscopy. Not only ground-state bands are studied, but there is also 
an increasing body of data for bands build on K isomers |2Ul I21| . Complementary information 
is obtained from decay spectroscopy of stopped ions. It provides detailed information about 7 
transitions from K isomers, a-decay transition energies, lifetimes, branching ratios for transitions 
to different levels in the daughter nuclei, and also information about electromagnetic transitions 
in the daughter nuclei. As in many cases the daughter nuclei are a emitters themselves, one 
experiment usually provides data for several nuclei [22l [23]. Finally, ion traps can be used to 
measure ground-state properties. The latter technique has been added only very recently to the 
arsenal of heavy-element research. The first pioneering experiments delivered the masses of a 
few nuclides around ^^^No |24| . Up to now, the quantum numbers of energy levels have been 
directly determined only for some heavy nuclei up to Z 100, cf. |151 [T8l [T9] and references 
therein. For heavier nuclei, the attribution of quantum numbers has to be based on arguments 
about selection rules for transitions, which sometimes leads to conflicting spin assignments or 
even level schemes when the experiments are performed with different set-ups and analysed by 
different groups. Again, the production cross section limits spectroscopic studies to the most 
favourable cases of projectile and target combinations. 

The experimental efforts are closely accompanied by theoretical studies. Mean-field models 
are the tools of choice for the description of the structure of transactinides. They come in two 
different forms, macroscopic-microscopic models on the one hand, and self- consistent mean-field 
models on the other hand. Macroscopic-microscopic models jHl [25| [26l [271 [28} [29] combine the 
" macroscopic" energy i?mac of a droplet of nuclear matter and a Strutinski shell correction -Emic 
that is calculated from the single-particle spectrum obtained from some parametrised single- 
particle potential. Both are calculated for a common shape parametrised by a set of coefficients 
{aj}, £'tot[{cei}] = -f'macil^i}] + -Emic [{««}] • There are many variants [30] that differ in the 
higher-order and deformation-dependent terms of the liquid drop and the form of the single- 
particle potential. A calculation consists in scanning the multi-dimensional energy surface for a 
given family of nuclear shapes for minima. By contrast, the entirely microscopic self-consistent 
mean-field models [31] . used in the form of HF-I-BCS or HFB, have as solution a set of single- 
particle wave functions {V'i} that minimises the binding energy i?tot[{^i}] from a given effective 
interaction and thereby provides local minima in the energy surface. Still, the surface can be 
mapped with the help of Lagrange multipliers. There are many variants that differ in the form 
of the effective interaction. The most widely used ones are the non-relativistic contact Skyrme 
IMl[Ml[Ml[35l[Ml[371[38l[39l[lQl[lIl[l2l[l3]and finite-range Gogny interactions [H [13 [Ml [ill [H] 
as well as relativistic mean-fields (RMF) l33[Ml[371[l9l[Ml[SIl[52l[53l[Ml[55l[56] that are 
also called covariant density functional theory. The macroscopic-microscopic method can be 
obtained as an approximation to self-consistent models [57]. In what follows, the concept of a 
shell correction will be used for the intuitive interpretation of the results from self-consistent 
mean-field calculations, even when it is not an ingredient of the model. 

Both methods have their specific advantages and disadvantages that are often complementary. 
The major technical advantages of the microscopic-macroscopic approaches are their numerical 
simplicity, that the parametrisations of the "macroscopic" and "microscopic" parts can be 
optimised separately, and that continuous multi-dimensional energy surfaces for fission or fusion 



can be easily constructed. The price to be paid is that this approach cannot describe the feedback 
of the single-particle wave functions onto the mean field, and that the convergence of the results 
in terms of higher-order shape parameters is difficult to control [28]. There is a multitude of 
parametrisations of the nuclear surface to be found in the literature [58], some of which are 
better adapted to compact shapes, whereas others are better suited for very elongated shapes 
and the transition to situations with two centres. The advantages of the self-consistent models 
are that due to their variational nature they automatically optimise all shape degrees of freedom 
not fixed by global symmetries and that they incorporate the feedback of the single-particle wave 
functions on the mean field. The latter is not only important for the description of polarisation 
effects driven by the Coulomb field, but also for the description of alignment effects in rotational 
bands and of the rearrangement of the nucleus when constructing K isomers. Extensions of self- 
consistent models can be used to calculate correlations such as quantum fluctuations in shape 
degrees of freedom [591 [60]. In addition, their time-dependent variants can be used for the study 
of certain aspects of the production of heavy elements [61]. Their main disadvantage is that 
the currently used effective interactions are not yet flexible enough to provide the same level of 
accuracy for the description of heavy nuclei as the macroscopic-microscopic models, and that 
the construction of continuous multi-dimensional deformation energy surfaces needed for fission 
studies with the help of constraints and Lagrange multipliers is technically challenging. 

The following sections will discuss three aspects of the structure of superheavy nuclei that are 
important for the extrapolation to the the island of stability based on the available spectroscopic 
information in the Z ~ 100 region. We will focus on self-consistent mean- field models. Section [2] 
presents predictions for single-particle spectra of spherical SHE, whereas Section [3] addresses 
deformed shell structure. Section jj] will explore the role of correlations beyond the mean field 
on the structure of SHE. Finally, Section [5] summarises the discussion. 

2. Spherical shell structure 

Shell structure of transactinide and superheavy nuclei is usually discussed in terms of the 
spectrum of eigenvalues of the single-particle Hamiltonian h 



in even-even nuclei. They reflect the mean fleld in the nucleus and provide a zeroth-order 
approximation to separation energies. It has to be stressed that there is no unique way of deflning 
a single-particle energy, and that most definitions are model-dependent. The shell structure of 
light nuclei is often discussed in terms of single-particle energies that are associated with the 
centroids of the strength function for nucleon pick-up and stripping reactions, cf. Ref. \Q2\ [63] 
and references therein. These cannot be expected to be equivalent to those from Eq. ([TJ. 

Even after more than forty years of experimental and theoretical efforts, the question where 
are the next closed spherical proton and neutron shells beyond ^"^^Pb is still open. There is no 
evidence that there are singly-magic spherical ones among the superheavy nuclei observed in 
experiment so far, and there is no consensus among theorists with regard to the location of the 
island of stability. Concerning predictions for shell structure, the main challenges for theory are: 

(i) The models have to be extrapolated far beyond the region where their parametrisations 
were adjusted. The reliability of the extrapolation of nuclear structure models has rarely 
been quantified so far. For mean-field models, flrst steps into this direction were taken just 
recently for the extrapolation to the neutron drip line at large proton-to- neutron ratio j64| . 

(ii) The average density of single-particle levels around the Fermi energy increases with particle 
number. As a consequence, the spacing of single-particle levels is much smaller than the 
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Figure 1. Proton single-particle spectra for 134 114 (left panel) and -'^20 (right panel) 
obtained with the folded- Yukawa single-particle potential used in macroscopic-microscopic 
models (FY), the Skyrme interactions SkP, SLy6, SkI4 and SkI3, and the NL-Z, NL-Z2 and 
NL3 parametrisations of the RMF. Data taken from Ref. j37j . 

one of lighter nuclei up to ^'^^Pb. Already a small relative shift of one single-particle level 
can open up or close a gap in the spectrum. 

(iii) The strong Coulomb fields try to push the protons as far apart from each other as possible, 
and thereby to polarise the nucleus. Corrections to the Coulomb, isovector and surface 
energies that play no evident role for lighter nuclei might become crucial for SHE. 

All of the above make predictions for the shell structure of superheavy nuclei in the island of 
stability a challenging task. The current situation can be summarised as follows: 

• All widely used parametrisations of the single-particle potentials for macroscopic-micro- 
scopic models (Nilsson, Folded- Yukawa, Woods-Saxon) predict Z = 114, = 184 as the 
dominant spherical shell closures in SHE [51 [271 I29j . However, modifying the 
surface thickness of the single-particle potential might close the Z = 114 gap and open up a 
gap at 126 instead [65l[66]- Note that the Woods-Saxon potential fitted by Chasman [151 [67] 
contains momentum-dependent terms that are not considered in the work of other groups. 

• Self-consistent models using Skyrme, Gogny, or RMF interactions almost always give 
different doubly-magic nuclei than the standard macroscopic-microscopic models. They 
predict either Z = 120 or 126 for the next proton shell closure and A^ = 172 or 184 for the 
next neutron shell closure |32[ [35l [37] . Not all combinations of these magic numbers are 
found, however. A notable exception is the SkI4 Skyrme parametrisation, which predicts 
Z = 114 and A^ = 184 like the macroscopic-microscopic models |35[ 137] . 

Typical results that illustrate the lack of agreement between the models for proton single-particle 
states are displayed in Fig. [H The comparison between f||ll4 and if2l20 illustrates also that 
the single-particle spectra from self-consistent mean-field models change rapidly with proton 
and neutron numbers in this mass region |35[ 137] . For example, the gap at Z = 120 is much 
larger when the proton number takes this value than for the Z = 114 isotope. A similar effect 
is found for the neutron shell closure at A^ = 172 [37]. As can bee seen in Fig. [H predictions 
for the gaps in the single-particle spectra rarely exceed 2 MeV. In fact, it has been pointed out 
that the actual gaps themselves are no longer the decisive criterion for magicity in SHE [39]. 
This is illustrated in Fig. [2] that shows the shell correction determined from fully self-consistent 
mean-field calculations with the Skyrme interaction SLy6 for the tin region and for superheavy 
nuclei. For nuclei in the Sn region, the combined shell corrections for protons and neutrons are 
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Figure 2. Total shell correction energy in 
MeV determined from fully self-consistent 
spherical calculations with the Skyrme 
SLy6 parametrisation for nuclei in the 
Sn region (bottom) and the transactinide 
region (top). Data taken from Ref. [39j . 



> 2 



@ 

^ / 2p3/2- 


lain k 'i'l^" 
126 1 




(§) 


^^_^;^;^=>-<.iV2- 



Figure 3. Proton (top) and neutron (bottom) 
single-particle levels of ^^^Sn (left) and 184I26 
(right) obtained with the Skyrme (SkP-SkI3) 
and RMF (NL-Z2, NL3) interactions as 
indicated. To illustrate the degeneracy of the 
levels, their thickness is drawn proportional to 
2j + 1. Data taken from Ref. [39] . 



localised at the magic numbers Z = 50, iV = 50 and N = 82. They take their largest values 
for the doubly-magic nuclei ^'^'^Sn and ^^'^Sn, and fall of quickly when going away from them. 
The situation is very different for superheavy nuclei: here, the shell correction takes large values 
in a large region that ranges from N ~ 164 to ~ 184 and Z ^ 114 to Z ^ 126, and that 
encompasses all possible predictions for doubly-magic superheavy spherical nuclei. This finding 
is fairly independent of the chosen effective interaction [39]. It can be understood when recalling 
that the spherical single-particle levels have a (2j -|-l)-fold degeneracy. The levels near the Fermi 
energy displayed in Fig. [T] range from j = 1/2 up to 15/2, and their degeneracy from 2 to 16. 
Its role can be better seen in the single-particle spectra displayed in Fig. O where the width of 
the levels is drawn proportional to their degeneracy. For spherical superheavy nuclei, the low-j 
orbits of protons and neutrons are bunched amidst a large number of high-j levels. Due to their 
large degeneracy, it are the latter that determine the average level density. The low-j levels can 
be shifted around opening up or closing the gaps between them without making a large difference 
for the level density and the size of the shell correction, even when the corresponding single- 
particle spectra look very different. The low-j are filled for proton numbers between 114 and 126 
and neutron numbers between 172 and 184, which is exactly the zone where the combined shell 
correction is large. The appearance of bunched low-j levels surrounded by high-j levels follows 
trivially from the sequence of single-particle levels in a finite potential well and the spin-orbit 
splitting in nuclei that is large enough to push an intruder into the major shell below. 

The appearance of the magic numbers Z = 120 and N = 172 in self-consistent models is 
the consequence of a polarisation of the density [371 1121 SH HHl [52] • The low-£ s, p, and d 
orbits that were found to be crucial for the extra binding of spherical superheavy nuclei in 
the discussion above have as a common feature that their radial wave functions are localised 
well inside the nucleus. In general, only \ow-i levels contribute to the density at small radii. 
Ignoring pairing and other correlations, ^72120 is the last nucleus for which they are not filled. 
By contrast, the occupied high-£ levels in the last major shell are all peaked on the the nuclear 
surface. Combining both, the net effect is a very pronounced central depression of both the 
proton and neutron densities, cf. Fig. [U that has been characterised as a "semi-bubble" density 
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Figure 4. Radial profile of the neutron 
(left), proton (right) and total density 
of the spherical configuration of f 72^20, 
obtained with the Skyrme interactions 
SkI3 and SkP, the RMF Lagrangian NL-Z 
and the folded- Yukawa (FY) single-particle 
potential widely used in macroscopic- 
microscopic calculations. Data taken from 
Ref. i37i. 




Figure 5. Axial deformation energy surface 
(top) and Nilsson diagrams of protons (middle) 
and neutrons (bottom) function of 

quadrupole deformation fi^-, Eq. ([2]), as 
obtained with the Skyrme interaction SLy6 for 
1^2^20. In the Nilsson diagrams, solid (dotted) 
lines denote levels of positive (negative) parity. 
The Fermi energies of protons and neutrons are 
indicated with dashed red lines. 



profile [HIHS]. In self-consistent models the density distribution feeds back onto the mean fields, 
in particular onto the radial profile of the spin-orbit potential |37j . For specific orbits such as 
the proton 3p and neutron 3d levels, the spin-orbit splitting is reduced to almost zero, cf. Fig. [H 
which opens up the Z = 120 and A'^ = 172 gaps. When filling the low-j orbits one goes back to a 
regular density distribution and the gaps at Z = 120 and N = 172 disappear [37] • Similar, but 
weaker, effects can be seen in lighter nuclei. For those it is usually just the non-occupation of one 
proton or neutron low-j level that causes a central depression of the density. The new feature of 
2^^120 is the much larger number of single-particle states and nucleons involved that makes the 
effect more pronounced. The size of this polarisation effect is correlated to the effective mass. It 
is the largest for small effective mass (NL-Z and SkI3 with mg/m ~ 0.6) and less pronounced for 
large effective mass (SkP, rriQ/m = 1). Interestingly, the density distribution from the folded- 
Yukawa (FY) single-particle potential frequently used in macroscopic-microscopic methods also 
exhibits a central depression, cf. Fig. [H However, in the absence of self-consistency it neither 
feeds back onto the shell structure, nor onto the macroscopic energy. 

Large gaps in the single-particle spectrum do not necessarily imply a spherical ground state. 
This is illustrated for ^y|l20 in Fig. [SJ The deformation energy surface is quite soft around the 
spherical point (/32 = 0), with an oblate minimum at a (32 value of —0.13. Its softness is again 
related to the presence of low-j levels around the Fermi energy for protons and neutrons. The 
splitting of their sub-states with deformation is rather weak and many of the levels on both 
sides of the Fermi surface repel each other such that the Z = 120 and A = 172 gaps persist far 
to the oblate side. The prediction whether ^^^120 has a spherical or a deformed ground state, 
however, is very sensitive to details of the effective interaction. The Gogny force [441 H5l HH] 
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Figure 6. Nilsson diagrams of proton (left) and neutron (right) single-particle levels of flpFniioo 
as obtained with the SLy4 Skyrme parametrisation [68] as a function of quadrupole deformation 
going from spherical shape to the prolate ground-state deformation. Labels on the left denote 
spherical j shells, labels on the right indicate jz and parity. Levels of positive (negative) parity 
are drawn as solid (dotted) lines. The colours indicate different mean values of j^. 

and most Skyrme interactions [A3\ I69j give a similar result as SLy6, whereas a few other Skyrme 
interactions [53] and most relativistic mean-field models give a spherical minimum [69j . 

There is a large number of non-equivalent definitions of quadrupole deformations to be found 
in the literature, which is a source of confusion and misunderstandings when comparing results of 
different groups. Throughout this article we will use the dimensionless quadrupole deformation 
f32 that is obtained from the mean value of the quadrupole moment taking out its trivial scaling 
with proton and neutron number 

where R = 1.2^4^/^ fm. Results from macroscopic-microscopic methods are usually shown in 
terms of the multipole expansion parameters of the nuclear surface, for which there exist several 
non-equivalent conventions and which are not equivalent to the mean values of the multipole 
moments [SH]. Results from self-consistent models have also sometimes been presented in terms 
of shape expansion parameters estimated from the multipole moments [62\ [31] . 

3. Deformed shell structure 

We have seen above that there exists a multitude of conflicting predictions for the shell structure 
of spherical superheavy nuclei that constitute the "island of stability" . As of today, the question 
about the location of the next shell closures cannot be answered by theory alone, and there is 
an ongoing effort to use the available spectroscopic data for well-deformed transactinides in the 
Z ~ 102 region. The basic idea can be illustrated by the Nilsson diagram shown in Fig. [HJ which 
covers deformations ranging from spherical shape to the prolate deformed ground state. At the 
ground-state deformation, the single-particle levels can be associated with the observed band 
heads of rotational bands in odd-^d nuclei. Measuring the energy differences between these band 
heads in the same nucleus and having a mean-field model that describes their relative positions, 
one can then estimate empirical spherical single-particle energies by tracing the Nilsson diagram 
back to spherical shape. The single-particle energies as they appear in the Nilsson diagram, 
however, are not observable. The differences between the energies of the various band heads 
measured differ from the single-particle energies in several respects. First, pairing correlations 
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Figure 7. Spectrum of eigenvalues of the mean- field (left, Eq. ([T])) and HFB Hamiltonian 
(middle, Eq. [3])) for protons and neutrons for the deformed ground state of ^^'^Fm as obtained 
with the SLy4 Skyrme parametrisation. Solid blue (dashed red) lines indicate eigenvalues of 
the HFB Hamiltonian corresponding to single-particle levels below (above) the Fermi energy. 
The right panel compares the spectrum of eigenvalues Eqp of the HFB Hamiltonian for neutrons 
obtained for a HFB vacuum constrained to Z = 100, N = 151 and the spectrum of the energies 
from self-consistently calculated blocked HFB states E'tot in ^^^Fm with experimental data for 
band heads taken from [TB]. All energies in the right panel are normalised to the ground state. 

introduce a correction to the relative energies. Second, the unpaired nucleon polarises the mean 
fields. In particular, due to the lifting of Kramers degeneracy in an odd-^ nucleus, there is 
no useful way of representing single-particle energies as defined in Eq. ([T|) for the purpose of 
analysing the shell structure. Third, each configuration in an odd-j4 nucleus has in general 
its own ground-state deformation. The effects of self-consistency when calculating these band 
heads are illustrated in Fig. [71 The left panel shows the single-particle spectra obtained for the 
ground state of ^^"^Fm, which correspond to the end of the Nilsson diagrams in Fig. O Pairing 
correlations modify their relative energies. These are incorporated in the spectrum of eigenvalues 
Eao of the HFB Hamiltonian 



often called quasi-particle energies and shown in the middle panel of Fig. [71 When neglecting off- 
diagonal matrix elements in the HFB equation dSj), which leads to the simpler HF-I-BCS scheme, 
the quasi-particle energy of a single-particle level represents its distance from the Fermi energy A 



and the size of its pairing gap A^^ through the relation E^p y (e^ — A)^ + A^^. For a realistic 
pairing interaction, the values of the gaps A^^ strongly depend on the single-particle level. In 
full HFB, the coupling of the quasi-particle states through non-diagonal matrix elements slightly 
modifies the spectrum. The right panel of Fig. [71 compares two different calculations of the band 
heads in ^^^Fm with data. One is a perturbative calculation of these states from the spectrum 
of quasi-particle energies E^p obtained from Eq. ([3]). The difference to the neutron spectrum 
in the middle panel of Fig. [71 is that the HFB vacuum is now calculated for particle numbers 
Z = 100, = 151. The readjustment of the Fermi energy for neutrons visibly rearranges the 
spectrum. In particular, the lowest quasi-particle state is now a 9/2~, in agreement with data. 
However, compared to experiment the spectrum of the Eqp is too compressed. Such perturbative 
calculation of low-lying states presents at best a crude approximation to the energies from full 
self-consistent calculations of the blocked HFB state for an odd-^ nucleus [TOl [7l], which is 
represented as Etot in the right panel. In this case, the energy of each configuration results 
from a separate calculation of the quasi-particle excitation performed on an optimised vacuum. 
There are many rearrangement effects, such as a reduction of pairing when blocking a state. 
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Figure 8. Single-particle spectra of neutrons (top panel) and protons (bottom) for 
the deformed ground states of = 152 isotones as a function of proton number (left) 
and of Hs {Z = 108) isotopes as a function of neutron number (right) calculated 
with the SLy4 parametrisation of Skyrme's interaction. 

and changes in the mean fields and deformations, which are also all different for each blocked 
configuration. The net effect is that the spectrum of the E'tot is markedly different: the level 
sequence changes, and the spectrum is more spread out. The calculated level sequence, however, 
does not agree with data. In particular, the 5/2"*" level should be lower and the 1/2"*" level higher. 
This can be achieved by closing the = 150 gap in the single-particle spectrum and making 
the = 152 gap larger. The same modification of the single-particle spectrum is suggested by 
the analysis of adjacent odd- A nuclei and by two-quasi-particle states observed as K isomers 
in even-even nuclei [18]. The presence of a deformed N = 150 gap not seen in experiment in 
connection with a too small = 152 gap is a common problem of all self-consistent mean-field 
models using Skyrme [SGl |l2] , Gogny [17j and relativistic mean- field interactions |36l EU [53] . 
Similar problems are found for the single-particle spectra of protons, where a gap should be 
at Z = 100, not Z = 98 and Z = 104 [22l [Ml Il2l 113 [SD [53]. The deformed gaps are fairly 
reproduced by the single-particle potentials used in macroscopic-microscopic calculations, which 
explains their superiority for describing details of the spectroscopic data in the A ~ 250 region. 

In any event, the single-particle spectra, in particular in the form of Nilsson diagrams, 
are a more intuitive tool for the interpretation of shell structure than the spectra from self- 
consistent blocked states. We will use them in what follows to illustrate two further difficulties 
of establishing the link between deformed and spherical shell structures. The first of these 
difficulties is that the analysis cannot be done globally using the band heads of nuclei throughout 
the region of heavy nuclei. Instead, it has to be done locally, using only data from one or 
perhaps a few adjacent nuclei. The reason is that there are substantial changes in the single- 
particle spectra of the deformed ground states when varying A^ and Z, cf. Fig. [TH The nucleons 
rearrange themselves for the maximum gain in binding energy that is possible for a given number 
of particles. In particular, some of the deformed shell closures such as Z = 116 and A^ = 162 
open up only for specific combinations of A^ and Z. The rearrangement of shells is accompanied 




Figure 9. Evolution of spherical proton 
shell structure as predicted by the Skyrme 
SLy4 parametrisation when going from heavy 
actinides to superheavy nuclei along the a- 
decay chain with N — Z = 50. For most nuclei 
shown, the spherical configuration corresponds 
to a local maximum of the energy surface, not a 
minimum. The spherical configuration of ^^'^Fm, 
i.e. P2 = in the left panel of Fig. [6l corresponds 
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by slight changes in all multipole moments and the density profile. The second difficulty is 
that also the spherical shell structure changes rapidly for heavy nuclei. This is illustrated by 
Fig. [9l Assuming for the moment that the effective interaction used would correctly describe 
the shell structure of nuclei around f5QFmioo, one would conjecture from the Nilsson diagram in 
Fig. [6] that there is a strong spherical Z = 126 shell closure. The ^ = 126 gap, however, shrinks 
substantially when going to heavier nuclei, which is the effect already discussed in Fig. [TJ In 
general, what is deduced about the magic numbers in the spherical single-particle spectrum in a 
non-magic nucleus might not be valid anymore for nuclei that actually have this magic number. 
Even when the SLy4 parametrisation does not describe all details of deformed shell structure in 
the Z ~ 100 region, the rearrangement effects are generic to all self-consistent mean-field models. 

The nucleon-number dependence of spherical shells as predicted by mean-field interactions is 
much stronger for SHE than for light elements, where the single-particle energies evolve much 
more slowly. In light nuclei, the observed strong quenching of various signatures of shell closures 
is only obtained when going beyond the self-consistent mean field \72 \ \73 \ IT U I75j. 

The deformed single-particle spectra are as sensitive to details of the effective interaction as 
are the spherical ones. Examples are given in Fig. [10] for Skyrme parametrisations with modified 
spin-orbit interactions. There are considerable changes relative to the Nilsson diagrams obtained 
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Figure 10. Parametrisation dependence of the Nilsson diagram of Fm for the example of 
three parametrisations of Skyrme's interaction. Left to right: SkI4 and SkI3 with two different 
non-standard choices of the isospin dependence of the spin-orbit field [76] and T44, which adds 
an explicit tensor force to a standard spin-orbit interaction |77j . 
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Figure 11. Kinematic (0*^^^ circles) and dy- 
namical (B(^\ diamonds) moments of iner- 
tia deduced from experiment (filled markers) 
and cranked HFB calculations with the SLy4 
Skyrme parametrisation (open markers). Cal- 
culated values taken from Ref. [l2], experimen- 
tal values from [JHl IH] (^sopm, 252,254]s^q) ^^^^ 
[791 [80] (256Rf). 



126 
122 
N 118 



correlation energy 








































1 4 


























































































































2 










































































































































1 















































































































































































136 110 111 118 152 156 160 164 168 172 176 180 184 



Neutron Number N 
Figure 12. Quadrupole correlation en- 
ergy in MeV from projection on angu- 
lar momentum J = and mixing of ax- 
ial configurations of different quadrupole 
deformation [78]. The correlation en- 
ergy is defined as the difference between 
the energy of the particle-number pro- 
jected mean-field ground state and the en- 
ergy of the particle-number and angular- 
momentum projected GCM ground state, 
cf. [73]. Calculations were done with the 
SLy4 Skyrme parametrisation. 



with SLy4, cf. Fig. [6l However, for none of these parametrisations the deformed = 152 gap 
opens up. One improvement is found for SkI4, which gives a deformed Z = 100 gaps like the 
macroscopic-microscopic models. This improves one- and two-quasi-particle spectra for proton 
excitations in the A ~ 250 region [7S]. However, it comes at the price of closing completely the 
N = 152 gap, which deteriorates the spectra for neutron excitations [75] . 

The proof of large ground-state deformations in this region of the nuclear chart has been 
delivered by the observation of highly collective rotational bands, first for ^^^No [ST], and 
then for some of its even-even and odd-74 neighbours [181 [T9] . The observation of these bands 
triggered numerous cranked HFB calculations, where the mean-field equations are solved with 
an auxiliary condition on the mean value of angular momentum [401 SD H2l H6l I51| . As an 
example. Fig. [11] compares calculations using the SLy4 Skyrme parametrisation with data for 
some even-even nuclei. The figure displays the kinematical (first) moment of inertia = J/lo 
and the dynamical (second) moment of inertia 0^^) = dJ/duj = Q^^^ + oj dQ^^^ /dio [181 ISS] as 
a function of the rotational frequency oj. The dynamical moment of inertia G^^^ is a sensitive 
probe for the internal rearrangement of the nucleons with increasing angular momentum J. The 
calculations agree reasonably well with the data, although there are disagreements in detail, for 
example the slightly different up-bend of G^^-* seen for ^^^No and ^^^No, which can be attributed 
to the deficiencies in the single-particle spectrum discussed above. 

4. Correlations beyond the mean field 

The mean- field approximation is best justified for heavy nuclei whose ground state corresponds 
to a deep minimum in a stiff deformation energy surface. This is indeed the case for many 
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Figure 13. Qa values of even-even nuclei 
from self-consistent mean-field and beyond- 
mean-field calculations |78j with the Skyrme 
interaction SLy4 compared with data where 
available |85j . Lines connect values for 
isotopic chains as indicated by the labels. 
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Figure 14. Dimensionless mass quadrupole 
moment /32 ([21) of the ground state of 
heavy nuclei as obtained with the Skyrme 
interaction SLy4. Positive values denote 
prolate, negative values oblate deformations, 
whereas /32 = indicates spherical shapes. 
Data taken from |78j . 



of the well-deformed mid-shell nuclei in the Z ~ 100, ~ 150 region that are the object of 
detailed spectroscopic studies. By contrast, most of the heaviest recently synthesised nuclides 
and their daughter nuclei with Z > 110, N >^ 164 are in a region of transitional nuclei |34j . 
which have soft deformation energy surfaces and in some cases also exhibit shape coexistence. 
For those, fluctuations in shape degrees of freedom or the mixing of coexisting shapes might have 
a significant impact on the ground state properties [721 [73l [TH [75l [56] . Both can be incorporated 
into the modelling in the framework of the generator coordinate method (GCM) |59t I60j . 

Figures [12] and [13] present results that consider pairing and quadrupole correlations beyond 
the mean field as described in Ref. [73]. The former are treated by projecting the mean- field state 
on good particle number. The nuclear shapes are restricted to axially symmetric and reflection- 
symmetric configurations. For these, the description of quadrupole correlations beyond the 
mean field simplifies from a five-dimensional calculation to the treatment of one Euler angle 
and the intrinsic deformation /32. The same Skyrme interaction SLy4 is used together with a 
density-dependent pairing interaction to generate the mean-field states and for the configuration 
mixing. A detailed analysis of the ground-state correlations obtained in such calculations for 
nuclei up to Z ~ 100 has been presented in Ref. [73j . These calculations were recently extended 
to superheavy nuclei [75]. The resulting beyond-mean- field ground-state quadrupole correlation 
energies are displayed in Fig. [T21 Their behaviour is correlated to the ground-state deformation, 
Fig-O similar to what is found for the transition from the rare-earth to the lead region discussed 
in Ref. [7H]. For the well-deformed nuclei below Z < 112 and .A^ < 166 the correlation energy is 
dominated by rotational energy from angular-momentum projection. Its value is on the order 
of 3 MeV and almost constant, and much smaller than the static mean- field deformation energy 
that is on the order of about 20 MeV. The correlation energy takes larger values up to about 
4 MeV for the transitional systems with 116 ^ Z < 124 and 168 ^ N < 184, and then drops 
rapidly to small values of about 1.2 MeV for the rigid spherical nuclei close to the double shell 
closure Z = 126, N = 184. For nuclei in the vicinity of spherical shell closures, the quadrupole 



correlation energy is on the same order or even larger than the static deformation energy. The 
variation of the correlation energy offers an explanation for the characteristic pattern of the 
measured two-nucleon separation energies across shell closures |74j in nuclei up to the Pb region, 
which is not described by pure mean- field calculations. Figure [13] illustrates its impact on Qa 
values of very heavy nuclei. The difference between the mean-field and beyond-mean-field results 
is evidently largest where the correlation energy displayed in Fig. [T2]varies quickly. This happens 
around the = 184 shell closure, where its variation from one nucleus to its daughter in an 
a-decay chain might be as large as 3 MeV, and to a lesser extend for nuclei just below Z = 126 
and throughout the transitional region down to Z values around 110. By contrast, for the well- 
deformed nuclei with Z ~ 100, the Qa values are not significantly affected as the correlation 
energy is almost constant. When drawing the Qa values for isotopic chains as a function of 
neutron number as in Fig. [T31 large gaps between the curves correspond to large deformed or 
spherical proton shell effects, whereas large downward jumps of the curves correspond to neutron 
shell effects. Both can be easily associated with the appearance of large gaps in the single- 
particle spectra displayed in Fig. [8l Again, the empirical shell effect at A'^ = 152 is not correctly 
described. This is a common deficiency of virtually all self-consistent mean-field models. On 
the mean-field level, very similar results are obtained also with other Skyrme interactions |36| . 
the Dls parametrisation of the Gogny force |47l HH] and with relativistic mean-fields |36l I50| . 
It is noteworthy that the Qa values from macroscopic-microscopic methods also underestimate 
the discontinuity at iV = 152, although they predict a sizable iV = 152 gap in the single-particle 
spectrum, cf. Ref. |301 183j . Curiously, the beyond- mean-field correlation energy has a visible 
impact on energy differences only around spherical shell closures, but not deformed ones. 

A similar method is the construction of a microscopic Bohr Hamiltonian through a series of 
approximations to projected GCM as described in Refs. (721 [75]. The study of selected a-decay 
chains using a five-dimensional collective Bohr Hamiltonian in the full (32-1 plane derived from 
a relativistic mean- field has been recently presented in Ref. |56j . 

For some of the nuclei discussed here, considering only axial and refiection-symmetric shapes 
might not be sufficient. For some of the transitional nuclei in the Z ~ 116, A'^ ~ 170 region, 
shallow triaxial minima have been predicted [3l]. Nuclei with neutron number slightly larger 
than 184 might be refiection asymmetric in their mean-field ground states [l8], which also hints 
at a possible octupole softness of adjacent nuclei. Likewise, based on the observation of low- 
lying 2~ levels at energies below two times the pairing gap it has been argued that nuclei around 

= 150 might be soft in non-axial octupole degrees of freedom |84j . 

An alternative for the description of correlations beyond the self-consistent mean field is the 
framework of particle- vibration coupling (PVC). There, correlations are added as a diagrammatic 
expansion in terms of single-particle degrees of freedom and RPA phonons. This methodology 
is currently restricted to spherical nuclei. A study addressing singly- and doubly-magic Z = 120 
isotopes in the relativistic mean-field framework |54^ finds that the single-particle strength 
becomes highly fragmented, consequence of the small gap in the unperturbed single-particle 
spectrum. The single-particle levels are shifted toward the Fermi energy, such that the gap 
in the single-particle spectrum becomes smaller. With 4-5 MeV, the PVC correlation energies 
of even-even Z = 120 nuclei reported in Ref. [55] are larger than the projected GCM values 
displayed in Fig. [T2j However, these values should not be directly compared, as the former 
include also correlation energy from pairing which is not contained in the latter. 

5. Conclusions 

To summarise the main points of our discussion 

• Because of the large density of single-particle levels, already small shifts of single-particle 
energies can qualitatively change the appearance of the single-particle spectra of SHE, which 
is not the case for lighter doubly-magic nuclei up to ^''^Pb. 



• For SHE, the importance of the spherical shell closures for magicity is different than in 
lighter nuclei. Because of the presence of highly-degenerate high-j levels and the smallness 
of the gaps in the single-particle spectrum, additional binding from quantal shell effects 
originates from the bunching of low-j orbits near the Fermi energy, not from the gaps. 

• It is not the magnitude of the shell effects that counts for the stabilisation of SHE, but how 
quickly it varies with deformation. One, sometimes even two, spherical shell closures do not 
automatically prevent superheavy nuclei from having a deformed mean-field ground state. 

• Self-consistency induces many local structural changes in SHE. 

• Their extrapolation to SHE amplifies differences between the models and parametrisations. 

• Almost all models agree that neutron shell closures enforce sphericity much more than 
proton shell closures do. This finding, however, has to be taken with a grain of salt. As 
pointed out in Refs. [73l[71], mean-field models tend to overestimate the size of neutron 
shell gaps for shell closures throughout the chart of nuclei up to = 126, whereas the 
proton shell closures are usually more satisfactorily described. 

• Quadrupole correlations beyond the mean field affect mass differences such as Qa values in 
the transitional region between the well-deformed actinides and the spherical shell closures. 
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